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ABSTRACT: This is a second paper describing the calculation of spectroscopy 

for orbitally excited states from lattice simulations of Quantum Chromodynamics. 

P_i| New features include higher statistics for P-wave systems and first results for the 
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spectroscopy of D-wave mesons and baryons, for relatively heavy quark masses. We 



parameterize the Coulomb gauge wave functions for P-wave and D-wave systems 
and compare them to those of their corresponding S-wave states. 



I. INTRODUCTION 

Recently we presented some preliminary results from an exploratory calculation of the 
masses of P-wave mesons and baryons from lattice Monte Carlo simulations of Quantum 
Chromodynamics in quenched approximation with Wilson fermions. This communication 
describes the final results from these simulations. We have increased the statistics of our 
P-wave study. There is a hint of some fine structure splitting in the charmonium system. 
We also present first results for D-wave spectroscopy of fairly heavy quark mesons and 
baryons. Finally, we show some of the properties of Coulomb gauge meson and baryon 
wave functions of P-wave and D-wave systems, and compare them to S-wave wave functions 
at the same quark masses. 

Only recently has lattice QCD spectroscopy begun to move beyond ground state 
hadrons. Some P-wave states' masses are regularly measured in staggered simulations 
because they are the odd parity partners of "ordinary" states: the a\ and p are examples 
of such pairs. In nonrelativistic QCD, Lepage and Thacker have computed the masses 
of xc an d XB states (without including spin effects). Few Wilson simulations have stud- 
ied P-wave states. The APE collaboration measured masses of some P-wave mesons in 
quenched simulations at 6/g = (3 = 5.7, but has had difficulty in continuing their program 
to higher (3. Recently, El-Khadra, Hockney, Kronfeld, and Mackenzie have presented a cal- 
culation of the IP-IS splitting in charmonium, which they use to fix the strong coupling 
constant. This calculation was done with a smaller lattice spacing than the one we report 
here, and with an improved action for the fermions. 

Calculations of the masses of orbital excitations in lattice simulations are difficult for 
three reasons: First, one needs to measure a correlation function with nonzero overlap onto 
the desired L sector and zero overlap on L = 0, otherwise the signal will be dominated 
at large t by the lighter L = states. Our methodology solves this problem. Second, the 
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signal is intrinsically noisy/ A diagonal correlator, C{t) = (0|r(t)r(0)|0), which falls off at 
large t like exp(— E\t) where E\ is the energy of the lightest state which the operator V 
can create from the vacuum, has fluctuations given by 

a2 = l((|r(t)r(o)| 2 )-c(t) 2 ). (i.i) 

Due to its first term, Op decays with a mass characteristic of the lightest particle |T| can 
make from the vacuum. If T is a meson operator (creating a qq pair) T will create a q q 
state, which will most likely couple to a tttt pair. Its correlator will fall like exp(—2m 7r t) . 
In the baryon sector |T| will make a q q state, and the lightest such state is three 
pions. Thus we expect a signal to noise ratio to be a falling function of t: u/Cn{t) ~ 
exp {mn ~ i7T-ir)t for mesons, and a/Cn(t) ~ exp {ran ~ 3/2m7r)t for baryons. This is 
a more serious problem for orbitally excited states than for S-wave states because their 
energy differences are larger. Finally, the baryon sector includes multiple states with the 
same quantum numbers, which will appear in the same correlators. For example, in the 
L = 1 [70] of SU(6)f the nonstrange sector includes one j = 5/2 and two j = 3/2 and 
j = 1/2 nucleons, and j = 3/2 and j = 1/2 A's. 



II. METHODOLOGY 



A. Construction of Orbitally Excited States 

We construct orbitally excited states by using interpolating fields which couple only 
to a specific angular momentum eigenstate, which are projected onto zero momentum and 
which are of large spatial extent to maximize overlap with the state. Our strategy is look 
at correlators of different operators at t = and t ^ 0. 

At the t 7^ end of the correlation function we use an operator which depends on the 
relative separation of the quarks, which is conventionally referred to as a "wave function. ' 
The wave function ipQ(r) of a hadron H in a gauge G is defined as 

V>G(O = ^(#k(i0^+r)|0> (2.1) 

X 

where q(x) and q(y) are quantum mechanical operators which create a quark and an 
antiquark at locations x and y. (We have suppressed Dirac and color indices.) Our 
correlation function is constructed from convolutions of quark and antiquark propagators 

G(x,y) 

C(r, t) = J>l*(i/1, y2)G q (m, 0; x, t)Gq(y 2 , 0;x + r, t)\0) (2.2) 

x 
where *&(yi, y<i) is the t — operator. At large t if the mass of the hadron is mjj , then 

C(r, t) ~ exp(— mjjt)ipQ(f) (2.3) 

and so by plotting C(f, t) as a function of r we can reconstruct the wave function up to 
an overall constant. We measure the mass of a state by convoluting C(f, t) with some test 
function which further projects out the desired state: 

C(t) = ^test(?)C(r,t). (2.4) 

r 
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At t — we take an operator which is separable in the coordinates of the quarks. For 
a meson we use 

V(xi,x 2 ) = faMfaM- (2.5) 

In order to couple to orbital excitations we take <p\ to be an S-wave and <p2 to be some 
orbitally excited state with angular momentum /, centered around some specified coordi- 
nate. This state is a linear superposition of a p = L = I orbital excitation and a state 
whose center of mass momentum is nonzero (this is the familiar "translation mode" of a 
shell model state). Convoluting quark propagators as in Eqn. (2.2) removes the p ^ 
state and gives us the wave function of the p = L = I state. 

B. Spin Considerations 

We did not construct a complete set of P-wave or D-wave mesons and baryons. Instead, 
we proceeded as follows (for P- waves): We worked in a basis in which 70 is diagonal. Our 
sources and sinks were chosen to couple only to the upper (large) components of the Dirac 
spinor. We constructed propagators for S-wave quarks with m s = ±1/2 and for P-wave 
quarks with mi = 1, m s = ±1/2. We can then completely construct the \jm) = 1 22) *P 2 
and \jm) = |11) P\ mesons, as well as the \jm) = |||) nucleon iV(5/2). We formed two 
other meson states with S = 1 and rrij = 1 and (mg = and -1) which are not angular 
momentum eigenstates; while we will label them as Pi and Po they couple to all three 
3 P states. (The " 3 Pi" state is (|| + IT)|H); the " 3 P " state is (||)|11).) If the ordering 
of states on the lattice were as in charmonium, the lightest state in the channel would have 
the smallest j for a given rrij and the wavefunctions would in fact project out the states 
which they label. In the baryon sector we constructed uud bound states with rrij = 3/2 
(TTI)lll) and -1/2 (111)111). 

D-wave spectroscopy is identical apart from the substitution of / = 2 for 1 = 1: we 
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completely construct the 3 L>3 and D2 meson states, and construct rrij = 2 and 1 states 
which overlap D2 and D\ states. We construct nucleons of rrij = 7/2, 3/2 and 1/2. 

This is an incomplete construction forced on us by computer memory limitations and a 
desire to keep the calculation simple. A dedicated simulation should do this properly, but 
will need many more quark propagators. Note also that while we are using nonrelativistic 
wave functions, they have the quantum numbers appropriate to the desired states and will 
couple to them, and not to S-wave states, regardless of whether the quarks are actually 
relativistic. 

C. Details of the Simulations 

We performed the simulations using the Connection Machine CM-2 at the Pittsburgh 
Supercomputing Center. Our data set consists of eighty lattices for P-wave studies, of 
which the last fifty were also analyzed for D-wave systems, computed in quenched approxi- 
mation at coupling (3 = 6, separated by a combination of microcanonical overrelaxed and 
Kennedy-Pendleton quasi-heat bath sweeps (100 passes of a pattern of four overrelaxed 
sweeps through the lattice and one heat bath sweep). The lattice size is 16 sites. We 
gauge fixed each of the lattice configurations to Coulomb gauge using an overrelaxation 
algorithm. We used Wilson fermions. We used a fast matrix inverter provided by C. 

1 *? 

Liu of Thinking Machines, Inc. to construct quark propagators. We computed P-wave 
spectroscopy at three hopping parameters corresponding to relatively heavy quark masses: 
k = 0.130, 0.1450, and 0.1520, and D-wave spectroscopy at ft = 0.1300 and 0.1450. The 
pseudoscalar mass in lattice units at these three k's is 1.43, 0.83, and 0.56. We used heavy 
quarks because the size of the wave function for the orbitally excited states are much larger 
than the size of an S-wave bound state and if the quark mass becomes too small, the wave 
function is squeezed by the simulation volume. 
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In Eqn. (2.4) we use a test function iptest — exp(— (r/rg) )r Y^ (O), in the relative 
coordinate. We take 4>\ in Eqn. (2.5) to be a Gaussian centered at the origin and </>2 to 
be exp(— (r/rg) )r Y, (O,). (As a technical note, when the width of the Gaussian becomes 
large compared to the size of the lattice one must worry about edge effects. We do that 
by replacing the spherical harmonic by a function which has the same symmetry but is 
periodic in the box (of size L): Y^ = (x + iy)/r^ (sm27rx/L + ism27ry/L)/r.) The width 
of the Gaussians used in \1/ and iptest was taken to be 2, 2, and 2v^2 lattice spacings for 
the three k's. 

We recorded wavefunctions for P25 P\i D3 an d D2 niesons and for JV(5/2) and 
N(7/2) baryons, at time slices 4, 5, and 6. All baryon wave functions pin the two quarks 
in a relative S-state to the same coordinate. We folded meson data onto one octant of the 
spatial lattice and baryon data onto one quadrant before storing it, and in addition kept 
data on one plane without folding. 



III. WAVE FUNCTIONS 



A. Pictures of Wave Functions 

Now we display some of the features of P-wave and D-wave wave functions and compare 
some of their simple observables to those of S-wave mesons of our earlier study. We begin 
our display of results for wave functions by showing some pictures of Coulomb gauge meson 
and baryon wave functions. 

In Figs. 1 and 2 we display plots of the real part of the wave function in the plane 
z = for a k = .152 iV(5/2) baryon and a k = .145 D% meson, respectively. The contours 
show the locations where ^ is a multiple of twenty per cent of its maximum value. The 
data for these graphs have not been spatially averaged. They show the characteristic 
dipole and quadrupole structure of the appropriate spherical harmonic. The fact that 
these distributions are not symmetric gives the reader an impression of the fluctuations in 
the data. They also give some idea of the extent to which the granularity of the lattice 
distorts the state, and the extent to which the state fits into the simulation volume. 

In order to further illustrate the extent to which a state fits in the simulation volume, 
we display in Fig. 3 a set of three dimensional contours of constant absolute value of 
the real part of the wave function. This figure, for the k = .145 D% meson, shows the 
characteristic four lobe quadrupole structure whose outer regions are compromised by the 
simulation volume. This is a problem for all the lighter mass states; all the k = .130 
mesons appear to fit reasonably well into the simulation volume. 

Finally, we extract the radial wave function f(r) itself, from 

f(r) = Re i)(r)/ReY}{yt). (3.1) 



The Coulomb gauge radial wave functions for k = .145 P2 an d -D3 mesons are shown in 
Figs. 4 and 5. They appear to show the characteristic linear and quadratic zeros of the 
wave function at the origin. 

B. Fitting Wave Function Parameters 

The goal of this section is to provide simple analytical parameterizations of wave func- 
tions which can be used for future studies of spectroscopy and to provide checks for calcu- 
lations of wave function properties performed directly on the data. 

Fitting the wave functions proved to be unexpectedly difficult because of the high 
correlations among wave functions at different separations. With only fifty or eighty lattices 
we had to fit a subset of the data since correlated fits require more lattices than fitted 
points. We elected to choose coordinates along axes where the spherical harmonic was 
unity (up to a sign); this gives us seven [z = 1 to 7) points to fit. We folded all directions 
related by reflections onto this axis with the necessary signs. Then we fit the radial wave 
function f(r) including a periodic or antiperiodic image term from the boundary (we fit 
ftest(r) = g{r) ± g(L - r)). 

After a certain amount of trial and error we chose to fit to 



g{r) = A(r + br ) exp(— cr) (3-2) 



for P- waves and 

g{r) = A{r + br ) exp(— cr) (3-3) 

for D-waves, plus image terms, and we believe that many other simple functional forms 
would work as well. The data cannot distinguish between these or more complicated 
functions, and if we try to force a fit, the Hessian matrix becomes singular. We fit the 
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baryon data to the same parameters (recall that we have pinned two quarks together, so 
there is only one relative coordinate left). 



The data are very correlated and the matrix of correlations very singular. It was not 
unusual to find correlation matrices whose conditioning number (ratio of largest to smallest 
eigenvalue) approached several thousand. The conditioning number was not stable: fitting 
half the lattices in a data set could cause the conditioning number to vary by a factor of 
two. We could also achieve a considerable variation in the conditioning number by varying 
elements in the correlation matrix by hand by a per cent or so. In contrast, the correlation 
matrices for propagators had conditioning numbers on the order of 50 and were quite stable 
under the same tests. Typically, in a 7 by 7 correlation matrix, only the largest 3 or 4 
eigenvalues remained reasonably stable as the number of lattices in the data set was varied. 
Therefore we adopted the following strategy for determining the parameters in f(r): we 
looked at uncorrelated fits, correlated fits to all parameters (very unstable), correlated fits 
in which the correlation matrix had its smallest three eigenvalues removed (via singular 
value decomposition) and correlated fits to a subset of the data (often r = 1,3,5,7). In the 
latter case one could not use consecutive points since either one or both of the correlation 
matrix or the Hessian matrix would become singular. The D-wave data was much more 
difficult to deal with than the P-wave data in this respect. 



The overall normalization of the wave function is not important for spectroscopy stud- 
ies. The parameters b and c for P-wave mesons are displayed in Fig. 6 and for D-wave 
mesons in Fig. 7. 
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C. Moments of Wave Functions 
The nth moment of the meson wave function is defined in terms of i^(r) as 

= MMiWW!. (3 . 4) 

J r^ar j(rj z 

The factor of 1/2 is included so that the second moment defined this way (when appro- 
priately weighted by quark charges) reduces to the second moment of the quark charge 
distribution defined through the form factor. 

We determined the first and second moments of our meson wave functions in two ways: 
first, we computed it directly from the data, by performing a single elimination jacknife 
analysis, and second, we computed it using the fitted form of wave functions. We consider 
the second method to be more reliable since in many cases the wave function is still large 
at the edge of the simulation volume. Only a fit which includes image effects can correctly 
reproduce the tail of the wave function. 

In all cases and for both jacknife and using the fitted radial wave function the values 
of two moments were independent of timeslice, although the uncertainty increases with 
increasing t. We display the first and second moments at t = 4 in Figs. 8 and 9 respectively. 
We see that while the two methods give rather similar results for the k = .1300 mesons, 
at smaller quark mass the discrepancy becomes pronounced. 

The P-wave wave functions are larger than the S-wave wave functions, and the D-wave 
ones are larger still. Notice that the diameter of the wave function in the simulation is 
four times (r) of Eqn. (3.4), so that the simulation volume we use would appear to be too 
small for D-wave systems made of lighter quarks. 
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IV. SPECTROSCOPY 

We extracted masses from our data by fitting the correlation function C(t) of Eqn. 
(5) in the standard way, and looked at "effective masses" (local slope of C(t)) and fits to 
a range t m \ n to t ma x = Tit/ 2 = 8. All data are fit including the effects of correlations at 
different times. 

As a general rule for selecting the best fit value to present in a figure or table, we use 
"fit histograms." A fit is represented by a rectangle centered on the best fit value for a 
mass /Li, with a width given by (twice) the uncertainty of the fit (i.e. \i ± A/i), and a height 
which is the confidence level of the fit (to emphasize good fits) times the number of degrees 
of freedom (to emphasize fits over big distance ranges) divided by the statistical error on 
the parameter (to emphasize fits with small errors). The fit with the greatest height is 
the one we quote. This was the method used to select the best mass in an earlier S-wave 
spectroscopy calculation. 

A. P-wave Spectroscopy 

Fig. 10 shows effective masses and fits to a range of t values for the meson and baryon 
data. 

All baryon masses at all k values appear to be consistent; there is little sign of a drift of 
the masses with choice of fitting range. There is no evidence of any fine structure splitting 
in any of the k values. One cannot say whether this is due to a small intrinsic splitting on 
the lattice, or whether all operators are merely coupling strongly to the j = 5/2 nucleon. 
In our extrapolations we will make the latter assumption. 

The k = .1300 mesons also have stable, consistent fits. The best fit values from 
histograms all begin at t = 2. There is a hint of the appearance of fine structure splitting 
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in the multiplet, as shown in Fig. 11. The splitting qualitatively resembles charmonium 
fine structure splitting, with the Pq state the lightest and the other states more nearly 
degenerate. However, uncertainties are so large that this probably should not be taken too 
seriously. These data seem to be limited by statistics. 

The k = .1450 data is noisier by about a factor of two. The P2 state is heavier than 
the P\. The Pq signal never stabilizes; while the fit from t m i n = 2 is satisfactory from 
the point of view of chi-squared (7.0 for 5 degrees of freedom) fits at increasing t m ^ n 's 
produce monotonically falling masses. Fit histograms are shown in Fig. 12. The P\ state 
is degenerate with the P2? but with large errors. 

Finally, the k = .1520 data share the same features as the k = .1450 data, with slightly 
larger uncertainties. The P-wave lattice masses are listed in Table I. 

B. D-wave spectroscopy 

The D-wave data is noisier than the P-wave data (as expected). Typical uncertainties 
for masses are about 0.08, four times the P-wave value. All masses appear to be asymptotic 
by tmin = 2 — 3 and all signals disappeared into the noise by t — 6. Because the data is so 
noisy, we could see no evidence for fine structure splitting in a multiplet. Fig. 13 shows 
effective masses and fits to a range of t values for the meson and baryon data. The masses 
are listed in Table II. 

C. Comparison to Experiment 

It is difficult to convert these lattice numbers into reliable quantities which can be 
compared to experiment. At (3 — 6.0 one is far from the scaling region. S-wave spectroscopy 
with Wilson fermions does not agree with experiment. The quark hopping parameters 
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we have used are very distant from the zero quark mass value. We will glance at two 
comparisons, but we have to say that with the quality of the signals they should probably 
not be taken very seriously as other than qualitative observations. 

The masses are shown in Fig. 14. First, if we extrapolate all masses linearly in k to 
k c = 0.1567, we find am( 3 P 2 ) = .93(3), ami 1 Pi) = .73(4), am(N(5/2)) = 1.07(4), and the 
D-wave meson and baryon are at 1.18(10) and 1.93(20) respectively. The proton mass at 
k c from our S-wave wave function study is am(N) = 0.55(1), so we have m( P 2 ) / m{N) = 
1.70(6), m( l Pi)/m(N) = 1.33(8), m{N(5/2))/m(N) = 1.94(8), and the ratios of the D- 
wave meson and baryon to the nucleon mass are 2.14(18) and 3.5(4). 

Experimental data corresponding to these states are m(a2)/m(N) = 1.40 or m(f2)/m(N) 
1.35, m{bl)/m{N) = 1.31, m(N (1675)) /m(N) = 1.78, m(p 3 (1690))/m(iV) = 1.80, and 
m(iV(2220))/m(iV) = 2.36. The P-wave masses look qualitatively correct, but the D-wave 
states (from which the extrapolation in k is enormous) are too high. 

As another comparison to experiment we can try to predict the mass of the D-wave 
states in charmonium. To do this, we must extrapolate in k from our k = .1300 data point 
to the charm mass. We also need a value for the lattice spacing a, which could vary by 
thirty per cent at this (3 value depending on how it is chosen. We determine k and a by 
taking lattice determinations of the P 2 and £1 states at k = 0.130, 0.145, and 0.152 and 
extrapolating their masses linearly in k. (We use the data of Ref. 16 for the k = .1300 
vector meson). We determine the lattice spacing by fitting the extrapolated masses to the 
V>(3095) and x(3555) masses. This gives a charm hopping parameter of k = .1224 and an 
inverse lattice spacing of 1/a = 1790 MeV. (Note that determining the lattice spacing from 
our proton mass data would give 1/a = 1710 MeV; from the rho meson, 1/a = 2264 MeV.) 
The extrapolated common D-meson mass is then 3.99(16) MeV, where the error is only 
from the extrapolation. The D\ cc state is at 3.77 GeV but its mass is influenced by the 
nearby DD threshold. Model calculations of D-wave states (some of which are narrow 
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since their decays to DD are forbidden) give masses of 3.81 to 3.84 GeV. At this value of 
the lattice spacing our P2 — Pq mass splitting at k = .1300 is 63 MeV; in charmonium 
the corresponding number is 145 MeV. 



V. CONCLUSIONS 

In retrospect, many aspects of the project could have been done better. We should 
have completely reconstructed the spin structure of all the different states. We should have 
recorded wave functions for all timeslices and not just for big Vs. Then we could have used 
the wave functions we determined in Sec. Ill as the input for i^test( r ) f° r spectroscopy 
calculations. 

It is clear that this program could be carried to arbitrarily high angular momentum 
states. To do so will probably require very high statistics, a more fine-grained lattice (since 
the lobe structure of the angular part of the wave function becomes more pronounced) 
and a larger simulation volume, since the size of the wave function grows with angular 
momentum. 

Note that the uncertainties on the P-wave masses fell by about a factor of two when 
the data set increased from 20 to 80 lattices. This suggests that the calculations of P-wave 
meson masses are almost certainly limited only by statistics. Meson P-wave spectroscopy 
still needs to improve its uncertainties by another factor of two to four before it can begin 
to make a serious comparison with experimental data, but we believe that this would be 
an easy thing for any large scale spectroscopy simulation to do. A reliable method for 
identifying specific baryon states remains to be demonstrated. 
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Table I. P-wave meson and baryon masses in lattice units. 



K 


state 


mass 


0.1300 


3 P2 


1.747(23) 


0.1300 


3 pi 


1.759(22) 


0.1300 


3 P0 


1.712(19) 


0.1300 


l Pl 


1.760(21) 


0.1450 


3 P2 


1.313(25) 


0.1450 


3 P1 


1.322(74) 


0.1450 


3 P0 


1.278(35) 


0.1450 


l Pl 


1.210(37) 


0.1520 


3 P2 


1.05(3) 


0.1520 


3 P1 


0.96(3) 


0.1520 


3 P0 


1.03(7) 


0.1520 


!P1 


0.88(4) 


0.1300 


baryon 


2.56(2) 


0.1450 


baryon 


1.79(4) 


0.1520 


baryon 


1.28(4) 



Table II.Dwave meson and baryon masses in lattice units. 



K 



0.1300 
0.1450 
0.1300 
0.1450 



state 



meson 
meson 
baryon 
baryon 



mass 



2.00(6) 
1.54(5) 
3.0(1) 
2-4(1) 
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FIGURE CAPTIONS 

1. Wave function of the k = 0.152 iV(5/2) baryon in the plane z — 0. The contours show 
interpolated lines of constant real ip in multiples of 0.2 times the maximum value of ip. 
The shading shows the value of the wave function (black is most negative, white most 
positive), interpolating from the original 16 lattice to a 64 grid. 

2. Wave function of the k = .145 D% meson in the plane z = 0. The contours and 
shading are parameterized as in Fig. 1. 

3. Surfaces of constant absolute amplitude of the k = .145 D% meson in three dimensions, 
in fractions of its maximum: (a) \ip\ = 0.2t(j m ax, (b) \ip\ = 0Aif; max , (c) \ip\ = 0.5ifj ma x, 
(d) IV'I = 0-Q^max- The "breaks" in the surfaces in (a) and (b) occur when the surfaces 
intersect the edge of the simulation volume. 

4. Radial wave function of the k = .145 P2 meson. Crosses show all points; squares are 
points along the x axis used to fit the wave function (fit is the line). 

5. Radial wave function of the k = .145 D% meson. Crosses show all points; squares are 
points used to fit the wave function (fit is the line). 

6. Fit parameters of P-wave wave functions (a) b parameter (b) c parameter. Points are 
labeled with crosses for P2 mesons, squares for P\ mesons, diamonds for iV(5/2) 
baryons. 

7. Fit parameters of D-wave wave functions (a) b parameter (b) c parameter. Points are 
labeled with crosses for D% mesons, squares for D2 mesons, diamonds for N(7/2) 
baryons. 

8. First moment of P-wave and D-wave meson wave functions (a) from a jacknife analysis 
and (b) from the fitted radial wave functions. Points are labeled with crosses for P25 
squares for P\, diamonds for -D3, octagons for D2 mesons, and fancy squares for 
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the pseudoscalar (data from Ref. 2). 

9. Second moment of P-wave and D-wave meson wave functions (a) from a jacknife anal- 
ysis and (b) from the fitted radial wave functions. Points are labeled as in Fig. 8. 

10. Spectroscopy of P-wave mesons and baryons. (a) meson effective masses, (b) meson 
fits to a range, (c) baryon effective masses and (d) baryon fits to a range. The three 
bands (in order of increasing mass) correspond to k = 0.152, 0.145, and 0.130. In each 
meson group, the Pi state is labelled by a cross, the Pq state by an octagon, the Pi 
state by a square, and P 2 state by a diamond. For baryons, the cross labels the P\/2 
state, the octagon the P 3 / 2 state, and the square the P 5 / 2 state. 

11. Fine structure splitting in the P-wave meson multiplets (a) k = 0.130, (b) k, = 0.145, 
(c) k = 0.152. 

12. "Fit histograms" of k = .145 P-wave mesons, (a) P 2 mesons, (b) Pi mesons, (c) Pq 
mesons, (d) Pi mesons. 

13. Spectroscopy of D-wave mesons and baryons. (a) meson effective masses, (b) meson 
fits to a range, (c) baryon effective masses and (d) baryon fits to a range. The two 
bands (in order of increasing mass) correspond to k = 0.145 and 0.130. In each meson 
group, the P> 2 state is labelled by a cross, the D\ state by an octagon, the L> 2 state 
by a square, and D3 state by a diamond. For baryons, the cross labels the P>\j2 state, 
the octagon the ^5/2 state, and the square the D 7 / 2 state. 

14. Masses of P-wave and D-wave hadrons as a function of hopping parameter. Diamonds 
label P 2 mesons, squares Pi, octagons Pq, and crosses, Pi. P-wave baryons are 
labeled with a fancy cross, D-wave mesons are labeled with a fancy square and D-wave 
baryons with a burst. 
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